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ABSTRACT 

The high stellar densities in the cores of globular clusters cause significant stellar interactions. These 
stellar interactions can produce close binary mass-transferring systems involving compact objects and 
their progeny, such as X-ray binaries and radio millisecond pulsars. Comparing the numbers of these 
systems and interaction rates in different clusters drives our understanding of how cluster parameters 
affect the production of close binaries. In this paper we estimate stellar encounter rates (T) for 124 
Galactic globular clusters based on observational data as opposed to the methods previously employed, 
which assumed "King-model" profiles for all clusters. By dcprojccting cluster surface brightness 
profiles to estimate luminosity density profiles, we treat "King-model" and "corc-coUapsed" clusters 
in the same way. In addition, we use Monte-Carlo simulations to investigate the effects of uncertainties 
in various observational parameters (distance, reddening, surface brightness) on F, producing the first 
catalog of GC stellar encounter rates with estimated errors. Comparing our results with published 
observations of likely products of stellar interactions (numbers of X-ray binaries, numbers of radio 
millisecond pulsars, and 7-ray luminosity) we find both clear correlations and some differences with 
published results. 

Subject headings: Globular clusters. X-ray binaries. Pulsars 



1. INTRODUCTION 

Soon after the discovery of bright X-ray binaries 
(XRBs) {Lx > lO'^* crg/s) in our Galaxy, it became 
apparent that they were overabundant (by a factor of 
'--^lOO per stellar mass) in globular clusters (GCs). This 
overabundance was att ributed to t he formation of XRBs 
by stellar interactions ()Clarklll975[ ). Models of how neu- 
tron star XRBs could be produced dynamically include 
tidal capture of a companion star by a neutron star 
(jPabian et al.lll975), coll isions of neutron stars with gi- 
ant stars ( Sutantvolll975l ) , and exchange of neutron stars 
into existing primordial binaries (Hills 197^). These in- 
teractions depend on bringing two stars, or a star and 
a binary, close together, and thus depend on the square 
of the stellar density. Gravitational focusing will bring 
stars closer together and is reduced by the stellar ve- 
locity dispersion, leading to a dependence of the stellar 
encounter rate (typically denoted T) on cluster proper- 
ties as r oc / p^/a, where p is the stellar density and cr 
is the velocity dispersion. 

Globular cluster stellar distributions have often been 
found to be accurately described by lowered, truncated 
Maxwellian potentials, known as King models (|Kingj 
[1961 [19661) . These models possess a core region of nearly 
constant and a rapid falling off of density outside the 
core. The majority of past work approximated the total 
r of a cluster by only considering the summed F within 
the core, assuming a constant density in the core; thus 
Fi oc p^r'^/a, where Vc is the physical radius of the cluster 
core. Additional approximations based on King model 
profiles have been used, particularly when a is not well- 
known for a clust er. In a King model p rofile, a cx p^'^rc, 
so F2 oc p^-^rl (jVerbunt fc Hutlll987| ). To date, even 
the most advanced calculations of F that have integrated 



p^ /<T have assumed the GCs follow a King model profile 
(e.g.,tooley ct al. 2003). 

These estimates have allowed comparison of the stellar 
interaction rates between different clusters in our galaxy, 
which showed that bright XRBs in Galactic globular clus- 
ters are indeed most concentrated in the highest-F clus- 
ters (|Verbunt fc Hut]ll987l: [Verbu nt 2003a). Ahhough it 
is more difficult to measure the surface brightness (SB) 
profiles of globular clusters in other galaxies, analysis of 
extragalactic globular cluster XRBs shows that they, too, 
tend to be conce ntrated in clusters that s how evidence of 
higher F values (I Jordan et al.l[200 4: Siv akoff et al.l 120071 : 
[Jordan et ani2007t iPeacock et al.ll2009l ). Evidence for a 
weaker-than-linear relation between F and the probabil- 
ity of hosting a bright X RB in other galaxie s (e.g., the 
nonlinear dependence of I Jordan et al.l [20041 can be ex- 
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plained by random err ors in the measurements of cl uster 
structural parameters. iMaccarone fc Peacoc3l2011| ). 

In our own Galaxy, however, we have accurate radial 
SB measurements of globular clusters, allowing precise 
estimates of F. The number of bright Galactic globu- 
lar cluster XRBs is still too small for precise tests of 
stellar encounter theories. However, recent X-ray, radio, 
and 7-ray observational advances provide la rge numbers 
of faint X-ray sources (jPoolcv & Hut''200l), radio mil- 
lisecond pulsars (MSPs, Ransom 2008), and integrated 
ray emission th at is presumed to arise from MSPs 
Abdo et al.ll2010[ ) . These results allow detailed compar- 
isons between F and the progeny of stellar encounters. 
X-ray binaries (both neutron star and white dwarf sys- 
tems) and millisecond pulsars (the descendants of X-ray 
binaries). 

However, current literature calculations of globular 
cluster stellar encounter rates only approximate the true 
density profile of the stellar cluster. The actual density 
profiles of many clusters do not exactly fit King models- 
there are 29 clusters in the Harris catalog (jHarrij 119961 
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2010 edition; hereafter HC) with designations of "core- 
cohapsed" , or possibly core-collapsed. Core-collapsed is 
an observational designation indicating that instead of 
showing a clear, flat central core, the radial SB profile 
of a cluster continues to increases towards its center. 
These observations are linked to theoretical models of 
a gravitat ional instability that lea ds to a rapidly shrink- 
ing core (jMevlan &: Heggij |1997[) . although the defini- 
tion of core-collapse used by theorists does not necessar- 
ily c oincide wit h the d efini tion used by observers (com - 
pare iHurlev fc Shara 2 Oil and iChatteriee et all |2012[ ). 
In addition to core-collapsed clusters, many clusters that 
have generally been considered to be well-fit by King 
models (e.g., NGC 6388) show radial SB gradients down 
to t heir centers, which are not predicted by King mod- 
els (jNovola fc Gebhardtl 120061 ). Calculations of stcl- 
lar encounter rate s usi ng different rnethod s (e.g. Fi, 
iHeinke et all 120031 : r9. iMaxweh et al.l 120121: integration 
of p^/tT of a King-model fit. iPoolev et aL 

MM) 

can get 

significantly different results, implying that the choice of 
method introduces a systematic uncertainty. This is a 
particular concern when considering how observationally 
core-collapsed clusters compare to other clusters, as none 
of the methods cited above use accurate descriptions of 
core-collapsed cluster properties (e.g., King-model fits to 
core-collapsed clusters simply assume a concentration pa- 
rameter, c, of 2.5, which overestimates the SB gradient 
outside the core). 

Moreover, previous calculations of stellar encounter 
rates have not, to our knowledge, quantified the uncer- 
tainties in their calculations. This makes it difficult to 
understand, when comparing F versus observations of 
close binaries, whether uncertainties in the input quan- 
tities, such as reddening, distance, or core radius (for Fi 
or F2), cause scatter in the correlations. 

Our goal in this paper is to rectify these two problems 
by calculating the 3-d radial luminosity density profile 
and integrating it to obtain an estimate of F. We then 
quantify the uncertainties in our calculations by Monte 
Carlo sampling from distributions of the observational 
inputs. Finally, we compare our results with some 
recent works to determine how our estimates affect the 
correlation of stellar encounter progeny with stellar 
encounter rates. Note that the goal of this paper is 
simply to perform an accurate computation of the 
simplest stellar encounter rate estimate, and its errors. 
We do not attempt here to include issues such as mass 
segregation, neutron star escape at birth, subsequent 
binary destruction, dynamical evolution of GCs, finite 
lifetimes, etc., which have been discussed in several 



1988 



IVerbuntI I2003bl: 
20081 ). as they do not 



works (e.g. iVe rbunt fc M evlanl 
iSmits et al.l 12006 : lyanova et al.l 
yet have simple, agreed-upon recipes that could be 
used to address these details. We will model these 
effects in an upcoming paper, where we will draw fur- 
ther conclusions about the dynamics of XRB production. 

2. DATA REDUCTION AND ANALYSIS 

To calculate F based on / (r) /a dr over several half- 
light radii, we need the luminosity density profile (as 
a function of radius), and velocity dispersions (ideally, 
also as a function of radius, but sec below), along with 
estimates of the distance modulus and extinction. 



2.1. Surface Brightness Profiles 

Our sample includes 124 Galactic GCs for which found 
published SB profiles. For 85 GC s we used the SB profiles 
compiled bv iTrager et al.l (|1995l hereafter T95). These 
datasets were obtained from various ground-based ob- 
servations, inastli;_Jiwn_theBeA Globular Cluster 
survey bv iDiorgovski fc KingI (jl986f ). T95 indicate the 
quality of the datapoints with a weight and their best 
d ata are labeled with w -eight = 1 . 

iNovola fc Gebhardtl ([20061 hereafter, NG06) provide 
SB profiles for 38 GC, some of which are also listed in 
HC. In these overlapping cases, we use the SB profiles 
provided by NG06 as they were constructed from Hub- 
ble Space Telescope (HST) observations, which are much 
higher resolution than ground-based data and were pro- 
cessed with attention to reducing the influence of the 
brightest (giant) stars. 

The quality of the observed SB data varies strongly 
from one GC to another (Fig. [ij. For all GCs except 
Terzan 5 (see details below), we used the Chebyshev 
polynomial fits provided in T95 or the spline fits pro- 
vided by NG06, instead of the raw photometric data. 
Given both the noise in the SB profile data and the strong 
dependence of our method on the derivatives of the SB 
profiles, we used the smoothed profiles throughout this 
paper. As we show in fJHl for GCs where the data is of 
high quality this choice has little effect on our calcula- 
tions. For GCs with poor quality data, the Chebyshev 
polynomial fits lead to a smoother luminosity density 
profile that should be more representative of the actual 
luminosity density profile. 

For three GCs (Palomar 10, Terzan 7, and Tonantz- 
intl a 2) the T95 SB profiles are uncalibrated. Follow- 
ing iMcLaughlin fc van der Marell (l2005l) . we calibrated 
these profiles by assuming that their central SB values are 
equal to the central SB values from the HC. For Terzan 
7, in addition to calibrating the data, we ignored the 
polynomial fit data for lograrcsec > 1.9 to avoid the non- 
physical increase of the fit SB profile with radius. Such 
a problem can be attributed to the lack of large-radius 
data points, and the high order of the Chebyshev polyno- 
mial fit. T95 also present two sets of data for NGC 2419. 
We choose the dataset which shows agreement with the 
central SB reported by HC. 

We estimated uncertainties on the SB profiles using the 
reported uncertainties in the photometric data. For the 
NG06 SB profiles, we used the maximum reported uncer- 
tainty in photometric data (requiring log rarcscc > 0). For 
the T95 SB profiles, we use d estimates of the photometric 
uncertainties calculated bv lMcLaughlin fc van der Marell 
(|2005h . 

As all the profiles were reported as a function of an- 
gular radius, we first calculated 1-D profiles as a func- 
tion of physical radius using the reported GC distances. 
To obtain 3-dimensional luminosity density profiles from 
the 1-dimensional observational SB profiles, we u s ed the 
non-parametric deprojection of iGebhardt et al.l (|1996f ) 
assuming spherical symmetry: 



p(r) 



1 r°° dp{R) dR 



dR 



(1) 
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Figure 1. Illustrations of SB profiles with different qualities. N fit is the fit provided by NG06. N photo are the photometric data points 
from NG06. T fit is the fit provided by T95. T photo are the photometric data points with weight=1.0 from T95. 



where fi{R) is the 1-D SB profile as a function of pro- 
jected radius (R) and p{r) is the luminosity density as a 
function of deprojected (spatial) radius (r). When cal- 
culating the luminosity density function, we first linearly 
interpolated the (T95 and NG06) fits to the SB profile 
to allow for a finer numerical integration. To integrate 
over the entire GC, we first set the central SB equal to 
the innermost data point (a very small extrapolation). 
We then set the integration upper limit to be the out- 
ermost available data point which is in all cases > 2.5 
half-light radii, checking to ensure that this truncation 
did not affect our final results. In some cases where the 
SBD decreases inside the core (e.g., due to noise or con- 
tribution of light from giant stars outside the core) , this 
integration yields a complex result. In all such cases, 
the imaginary component is less than 10~^ the size of 
the real component. By ignoring this small imaginary 
component, we can reliably calculate the radial density 
distribution. Fig. [2] shows the result of interpolation and 
deprojection for NGC 104, one of the most well-studied 
clusters. 




10° 10^ 10^ 10^ 

Radius {arcsec] 



Figure 2. Interpolated surface brightness profile (top, interpo- 
lated from values given by NG06) and deprojected luminosity den- 
sity profile (bottom) for NGC 104. 
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Figure 3. Comparison of tiie modified central velocity dispersions 

from G02 (o"Modificd Gnedin — '^Gnedin 

/1. 57) to those from HC for 
62 GCs in common (values in km/s). The line represents auc = 
fMG- 
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Figure 4. Velocity dispersion profile for NGC 104. The solid line 
is the projected profile, the dashed line is the deprojected profile 
obtained from the sum, Eq.|4l and the dotted line is the deprojected 
profile obtained from the integration, Eq. [3] Note that the core 
radius for NGC 104 is Ri 21". 



2.2. Velocity Dispersion 

We have full velocity dispersion (ct) profiles for only 14 
clusters (see Table [3] for these sources). For the remain- 
ing clusters, we only consider the central <j value. Since 
a falls off much more slowly than SB with radius in the 
cluster, using the central a value for all radii produces 
very small changes in the inferred V (see Our pri- 
mary source for central values of a and their errors was 
HC, which compiles central velocity dispersion measure- 
ments for 62 GCs (gnc)- For other GCs, we referred to 
theoretical estimates by IGnedin et al.l (|2002L hereafter 
G02). For GCs where HC reports velocity dispersion, 
the G02 values are 1.57 times larger on average. So for 
the cases where HC does not report velocity dispersion, 
we used modified values from G02, (ctmg = crG02/l-57). 
Fig. [3] shows our comparison between the a values from 
G02 (modified) and HC for the 62 clusters in common. 

For GCs where HC reports velocity dispersion, we used 
estimations he provides for uncertainty in the velocity 
dispersion. For the rest of our sample, we used the aver- 
age fractional discrepancy between ctmg a-nd cthc for the 
62 GCs they both report, as our uncertainty 5: 




0.32 



(2) 



For the 14 GCs where we had detailed velocity disper- 
sion profiles, we could compare the effects of assuming 
a constant velocity dispersion instead of using the true 
velocity dispersion profile. For these clusters, we depro- 
jected the ID velocity dispersion profile to a 3D profile 
making the assumption of spherical symmetry. We used 
the non-parametric integration for deprojection: 



p{r)a{r) 



1 



d{^l{R)(7.p{R)) dR 



dR 



Vr^ 



(3) 



where ap{R) is the projected ID profile and a-{r) is the 
deprojected 3D profile. Since the velocity dispersion data 



had not been previously smoothed, we applied a third- 
order interpolation prior to deprojecting the velocity dis- 
persion. We truncated the integration at the outermost 
data point. This method produces a drop to zero at 
the outer radii, due to our choice of integration limits 
(choosing the outermost data point instead of infinity). 
To check the overall validity of the first method, we used 
a second method of deprojection assuming a discrete sum 
of shells, where we set CTp and a to be equal in the out- 
ermost layer of the GC (we omit the factor of -s/S in 
converting from 1-D to 3-D velocities, as it will be iden- 
tical in all clusters, assuming isotropic orbits). By these 
assumptions, we calculate a discrete sum for the projec- 
tion: 

aARn) - ^^^^ (4) 

where i starts from the outermost radius and goes to- 
wards the center. Starting from the outermost layer, we 
found values for the deprojected a at different points and 
interpolated them as a function of r. To compute F for 
these 14 GCs, we used the deprojected profile obtained 
from the latter method. In Fig. |4]we present a compar- 
ison of the projected dp profile and the deprojected a 
obtained from both methods for NGC 104. In fJHwc dis- 
cuss the effects on F of using a full deprojected a profile 
versus assuming a constant tr throughout the cluster. 

2.3. Distance Modulus and Extinction 

To estimate luminosity density as a function of phys- 
ical radius for GCs, we need to calculate the physical 
radius using distance and angular radius. To calculate 
distance and estimate uncertainties on it, we used val- 
ues for the apparent distance modulus D[m — M) and 
foreground reddening E{B — V) from HC. Based on the 
different claimed measurements in the literature for a few 
GCs, we assumed an uncertainty of 0.1 magnitude in dis- 
tance modulus for all GCs (except Terzan 5, see below). 
Following HC, we generally assumed a 10 % uncertainty 
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Figure 5. Available SB data for Terzan 5. T fit: fit provided 
by T95. T photo: photometric data points with weight=1.0 from 
T95. L photo: Photometric data from LIO. In comparison LIO 
shows a higher SB value for the core. 



for the reddening, E{B — V), imposing a minimum un- 
certainty of 0.01 magnitude for any cluster. We used 
Ay = RvE{B — V) with Rv = 3.1 to obtainthe extinc- 
tion. Since Ry is not the same for all p arts of the sky 
([Hendricks et al.ll2012l: INataf et al.ll2012l) . we assumed a 
further uncertainty of 10 % in Ry. For 3 GCs (AM 1, 
NGC 5466, and NGC 7492) HC reports E{B ~V) ^ 0.0, 
in these cases we used alternative source s to improve 
these estimates. For AM 1 we chose 0.02 (iDotter et al.l 
[2008I) . for NGC 5 466. 0.02 dSchlegel et al.l ll998D. and for 
NGC 7492, 0.04 (ISchleeel et al.lll998l) . 

2.4. Special case of Terzan 5 

Terzan 5 is a highly extincted GC near the Galac- 
tic core that contains > 50 XRBs (iHeinke et al.ll2006D 
and > 33 millisecond radio pulsars (jRansom et al.ll2005 ) . 
This large population of sources makes it an ideal GC for 
more detailed analysis. Although SB profiles are avail- 
able i n T95, we note that h igher quality data was avail- 
able in lLanzoni et al.l ()2010l hereafter LIO), derived using 
7J5Tobservations (ACS - F606W). However, LIO did not 
provide clear fit parameters. As a result, we use their 
photometric data to derive SB (Fig. [5]) . We assume an 
uncertainty of 0.2 magni tudes for the S B pro file, as re- 
ported by LIO. Recentlv iMassari et al.l ()2012D presented 
a high resolution reddening map of Terzan 5. From their 
map, we find E{B — V) ~ 2.61 for the core of Terzan 5 
and used their estimate of Ry = 2.83 to obtain our Ay 
estimate. For its distance modulus we used the value 
of (m — M)„=21.27 from H C, which with our A y gives 
the same (m— M)o=13.87 as lValenti et al.l ()2007[) derive. 
However, due to the uncertainty in measuring (m — M)o 
in this highly reddened case, we assumed a conservative 
uncertainty of 0.2 for this quantity. 



3. STELLAR ENCOUNTER RATE, T 

To calculate F, we numerically integrated ap{r)'^ / (T{r) 
using the luminosity density and velocity dispersion pro- 
files derived above, where a is an arbitrary constant that 



was set by requiring the F value for NGC 104 be equal to 
1000. To ensure that the first-order interpolation of the 
fits to the SB profile were appropriate, we recalculated 
F using both second-order and third-order interpolation. 
This led to no significant differences in the final results 
(< 0.1 % change). 

To estimate the uncertainty in F, we performed Monte- 
Carlo simulations of the F calculation with different in- 
puts. Our principal code is written in Mathematic^ 
and the average number of iterations for each GC was 
« 400. We assumed gaussian distributions for the in- 
put parameters (distance modulus, reddening, Ry, SB 
profile amplitude, velocity dispersion) with the reported 
values as the mean value, and the reported uncertainties 
as the standard deviations of the distributions. We used 
these distributions with caution, modifying them when 
they were unphysical. For low values of extinction, the 
gaussian distributions include negative values. For the 
velocity dispersion, values very near to zero also produce 
unphysical results (since velocity dispersion is in the de- 
nominator in F ~ J p/a). So we did not run simulation 
for those values. 

In the case of extinction, we required positive values, 
and in the case of velocity dispersion, we required that 
the simulated velocity dispersion was within two stan- 
dard deviations (eqU)) of the measured velocity disper- 
sion. For 2 GCs, NGC 7492 and NGC 5946, the reported 
uncertainties from HC on a are more than 50 %, so for 
these two, we truncated the a distribution at one stan- 
dard deviation instead. 

When the photometric data was of high quality, we 
found that integrating this data directly gave similar re- 
sults as integrating the fitted Chebyshev polynomials. 
The differences in the final results were typically < 1 % 
(e.g. NGC 104). In Table [T] we provide a comparison be- 
tween F calculated based on the photometric data, and 
based on the Chebyshev fit for some of the GCs where 
data were available from NG06. In the few cases with a 
large difference between the two values (e.g., NGC 5897, 
NGC 6205 & NGC 6254), the observational data did not 
extend out to the outer portions of the GC. In these 
cases, by truncating the Chebyshev fit profile to the out- 
ermost point of the photometric data, we greatly reduce 
the difference in results; for NGC 5897 it drops to 26.7 % 
and for NGC 6254 to 13.2 %. 

4. RESULTS 

The final F values we report (Table [2]) are calculated 
based on the default values for quantities described in 
^ For most clusters, the F calculated from the default 
values lies within 5% of the median of the histogram of 
F values produced in our simulations (Generally the dis- 
crepancy between default value and median of the distri- 
bution is caused by truncation of the input parameters 
distribution described in ^ . Uncertainties in F for each 
source are calculated based on the histograms of F val- 
ues produced from our Monte-Carlo simulations (Fig. [5]). 
We identify the 1-cr upper bound by increasing F from 
the median of the distribution upwards until we include 
an additional 34% of the simulations, and similarly iden- 
tify the l-CT lower bound. (Note that the F probability 
distribution is not necessarily a Gaussian.) 

^ http://www.wolfram.com 
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4.3 




172.4 


194.4 


11.3 




449.0 


458.0 


1.9 




207.1 


191.1 


8.3 




1046.4 


984.3 


6.3 


NGC 5897 


0.2845 


0.850 


66.5* 


NGC 5904 


152.42 


164.1 


7.1 


NGC 6093 


568.24 


531.6 


6.9 


NGC 6205 


48.475 


68.91 


29.6* 


NGC 6254 


13.656 


31.37 


56.5* 


NGC 6266 


1827.1 


1666.5 


9.6 


NGC 6284 


670.77 


665.54 


0.8 



Table 1 

Comparison between F calculated based on photometric data and 
based on the Chebyshev fits to the SB profiles (from NG06). 
Incompleteness in the photometric data appears to explain the 
cases with a large difference between the two F values (marked by 
a *). In these cases, truncating the fit to the region where 
photometric data is available reduces the difference in results. 




100 700 aOO 900 1000 1100 1200 1300 1400 1500 
Normalized Stellar encounter rate 

Figure 6. Results of Monte-Carlo simulations for NGC 104, show- 
ing the number of trials giving each value for F. Due to our choice 
of normalization, the histogram is forced to be centered on 1000. 



Name 


F 


Lower bound 


Upper bound 


Terzan 5 


6800 


3780 


7840 


NGC 7078 


4510 


3520 


5870 


NGC 6715 


2520 


2250 


2750 


Terzan 6 


2470 


753 


7540 


NGC 6441 


2300 


1660 


3270 


NGC 6266 


1670 


1100 


2380 


NGC 1851 


1530 


1340 


1730 


NGC 6440 


1400 


923 


2030 


NGC 6624 


1150 


972 


1260 


NGC 6681 


1040 


848 


1310 


NGC 104 


1000 


866 


1150 



Table 2 

The 11 GCs with the highest F values, providing their F values 
(normalized to give NGC 104's F=1000), and 1-cr bounds. The 
complete table is available in the electronic edition of the journal. 
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Figure 7. Comparing Fi = p^r^cic^ and F = a^^ f p^{r)<fir. 
Core-Collapsed GCs are denoted in red. 



name 


Difference (%) 


Ref. 


NGC 104 


1.13 


(1) 


NGC 288 


2.86 


(1) 


NGC 362 


2.63 


(1) 


NGC 2419 


2.94 


(1) 


NGC 3201 


5.00 


(1) 


NGC 5024 


0.23 


(2) 


NGC 5139 


1.32 


(1) 


NGC 6121 


0.08 


(1) 


NGC 6218 


14.7 


(1) 


NGC 6254 


0.15 


(1) 


NGC 6341 


2.56 


(1) 


NGC 6656 


4.05 


(1) 


NGC 6809 


7.89 


(1) 


NGC 7078 


2.84 


(3) 


Table 3 



Difference between F calculated based on a constant velocity 

dispersion and the measured velocity disp ersion profile. 

References for velocity dispersion profiles - (1): IZocchi e t al.ll2012l 
(using their Kin g mod el fits to the profiles), ( 2) ]3oUima'etan 
[20ig. (3') IMurphv et al.|[201ll 



We also investigated the effects of assuming a constant 
velocity dispersion profile by comparing the F computed 
based on a constant a profile versus the actual measured 
(and deprojected) a profile for 14 GCs. For the purposes 
of this comparison alone, we used the central velocity 
dispersion values reported by these profiles as the value 
for the constant velocity dispersion calculations (instead 
of the values from HC or G02). For deprojecting the 
observed velocity dispersion profiles we used the method 
of sums described in Section 12.21 As shown in Table 
[31 the difference between the results is always less than 
15%, and usually less than 5%. For the final values, 
for consistency, we used a constant a for all GC for the 
calculations presented in Table [2] and El 

To have a complete set of calculations, we also included 
F calculations and uncertainty estimations based on the 
simplified equations {plr^<Tc^ pI'^'^'I) ^ot 143 GCs 
(19 GCs in addition to the main sample) using HC (Ta- 
ble [5]) . To do this, we used central surface brightness (in 
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magnitude per arcsec^), extinction, distance modulus, 
estimated core radius, Tc and concentration parameter, 
c, to calculate centra l luminosity densi ty, pc- Following 
the prescription from iDjorgovskil ()1993f ) : 



-[^g0.4(26.362-^i) 

P = 

P rc 

_ in-0. 603x10^^+0.302 



(5) 



where p = iQ-^> bOixiD +0.302 jg parsec. 

For velocity dispersion we used the central values that 
we aggregated from the literature in ^2.21 Similar to the 
method described in S|3]we did Monte-Carlo simulations 
to estimate their effects on F. We assumed uncertainties 
on extinction, distance modulus and surface brightness 
as before. We also assumed an uncertainty of 5% for 
the core radius. Since the concentration parameter c has 
little effect, we did not include any error on c. Compar- 
ing these simplified values of F to our main results, the 
differences are relatively small for many GCs (Fig. [7]). 
Although the value of F for Terzan 10 calculated by the 
simplified method (F2) is extremely high, we found it 
to be untrustworthy. While HC reports the core radius 
of Terzan 10 is ^ 0.9', inspection of a 2-MASS J-band 
image from the Infrared Science Archival suggests it is 
< 0.2'. 

5. APPLICATIONS 
5.1. X-ray sources 

A significant difference between our results and pre- 
vious works comes in the case of c ore-collapsed clusters. 
For instance. iMaxwell et al.l (|2012[ ) derives similar values 
for F2, with differences principally arising in the core- 
collapsed clusters (Fig. |8]). 

Comparing our values for F to results from 
iPoolev et al.l (|2003| ) (which calculate F by integration 
over the half- mass radius assuming king models) and 
lFregeaul([200l . our calculations show that, at about the 
same values of F, core-collapsed GCs have lower numbers 
of X-ray sources compared to t ypical GC s (Fig. [H) . This 
is in c ontrast with the results of lFregeciul (j2008). Fregeaul 
(|2008f ) suggested that, contrary to previous thinking, 
most globular clusters are currently still in their "early" 
contraction phase, and that only those clusters obser- 
vationally defined as "core-collapsed" have reached the 
binary-burning phase. These clusters would then need 
to be currently "burning" binaries to support themselves 
at their current core radius. The initial impetus for this 
suggestion was the apparent excess of X-ray sources in 
three "core-collapsed" clusters, NGC 6397, M30, and 
Terzan 1, compared to other GCs with similar values of 
F. This would be explained if X-ray binaries were created 
a few Gyrs ago, at a time when non-corc-collapsed clus- 
ters were substantially larger and less dense, but core- 
collapsed clusters presumably were at their current size. 
Thus, the F relevant for producing the current X-ray 
sources in non-core-collapsed clusters would be smaller 
than the currently observed F, as those clusters will have 
contracted and become denser. Our calculations remove 
the evidence for NGC 6397 and M30 having higher-than 
average X-ray source numbers for their F values. Instead 
our results suggest that core-collapsed clusters underpro- 
duce X-ray binaries. 

^ |http: //Irsa. ipac . caltech. edu/| 



NGC6266 
NGC6388 \ 



NGC6626 
NGC6093^ 
NGC6752 7^ 
^^'^ 



NGC104 



NGC5904 
NGC5139 NGC7099 
NGC639r 



NGC6121 



r This work 

Figure 8. Comparison of F estimates by Maxwell et al. 2012 ver- 
sus values calculated in this work (using a different normalization) . 
Corc-coUapsed clusters have errors shown in red, and show many 
of the largest differences. Note also that NGC 6388 has a lower F 
in our calculations. We choose our F normalization to give average 
values similar to those of Majcwell et al. 2012. 



One cluster that may not fit with this picture is Terzan 
1. This is a GC that appears to be core collapsed, but its 
structural parameters are poorly determined at present. 
However, its position near the Galactic centre suggests 
an altern ative scenario , that it may have been tidally 
stripped (|Cackett et al.l 120061 ). 

Fig. [9] indicates that, altho ugh the a sserti on about 
globular cluster evolution by iFregeaul (|2008| ) may or 
may not be true, the numbers of X-ray sources above 
4 X 10'^'^ erg do not provide evidence for this assertion. 
Other evidence, perhaps from comparing detailed Monte 
Carlo models of gravitational in teractions between star s 
with observed quantities (e.g., iChatteriee et al.l |2012|) . 
may illuminate this question. On the other hand, the 
X-ray sources in core-collapsed clusters will ex perience 
substantial binary destruction (jVerbuntI l2003b( ). which 
may explain their r a ther different l umino s ity functions 
(IPoolev et al.l 120021: iHeinke et ail l2003t iStacev et al] 

5.2. Numbers of radio MSPs 

Large numbers of radio MSPs have been detected in 
several GCs, with the largest numbers in very high-F 
clusters (jCamilo fc Rasiol [20051 ). Several works have at- 
tempted to compare the numbers of MSPs in different 
clusters, accounting for the detection limits of the surveys 



to MSP numbers ( 


Johnston et al.l 119921: iHessels et al. 


2007 


; RansomI 12008 


: Hui et al.l 120101: iLvnch & Ransom 


2011 


; Bagchi et al.l 


201 If). These analyses must esti- 



mate the radio luminosity functions of cluster MSPs 
and the sensitivity of different surveys (involving com- 
plex estimates of pulsar detectabili ty). Perhaps t h e mos t 
sophisticated of these is that of iBagchi et al.l (|2011| ). 
which incorpor ates information from diffuse radio flux 
measurements (|Fruchter fc Gossll2000HMcConnell et al l 
20041) a nd summed 7-ray emission (e.g. lAbdo et al.ll2'oTor 
Hui et al..,2011i : see also below). 
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! 10' 



NGC5139 
NGC6397_ 



NGC7099 
T NGC6752 
NGC5904 ^|-|^ — I 



Stellar encounter rate r 

Figure 9. Numbers of XRBs (from Poolcy et al. 2003, Lugger et 
al. 2007), excluding the GCs with only lower limits determined), 
compared to our F estimates, with appropriate error bars. Core- 
collapsed clusters are in red, and show a tendency towards fewer 
XRBs for their T than other clusters. We choose a T normalization 
assuming r=20 for NGC 7099. 



10" 10' 10^ 10^ 

stellar encounter rate 

Figure 10. Number of recycled pulsars within a GC (from 
Bagchi et al. 2011) compared to our calculated T values. Core- 
collapsed clusters are shown in red. A correlation is clearly seen, 
and tabulated in Table 2] The normalization of F is chosen to 
bo similar to the choice of normalization in Bagchi et al. 2011 
(rNGC 6266 = 100). 



IBagchi et al.l (j2011[) calculate the most likely numbers 
of MSPs in 10 globular clusters, based on their simula- 
tions of the detectability of MSPs in these clusters, and 
from the observations discussed above. They make the 
striking claim that there is no compelling evidence for 
any direct relationship between any GC parameter and 
the number of MSPs per cluster; in particular, they claim 
that th ere is no correlation between F and the number of 
MSPs. IBagchi et al.l (I20T1I) use Pearson, Spearman, and 
Kendall statistical correlation tests, and report the rel- 
evant coefficients and null-hypothesis probabilities. We 
note that the null-hypothesis probabilities for the Spear- 
man and Kendall tests for correlation between their cal- 
culated r and the numbers of MSPs are 0.02 and 0.01, 
rather less than the typical 0.05 criterion for significance. 
However, the Pearson test's null-hypothesis probability 
is only 0.07, which does not provide clear evidence of 
correlation. 

Here we assume that their calculations of the num- 
bers of MSPs are correct, and recalculate these correla- 
tions usmg_ou£MwJ]values. We use model 1 (FK06) 
from lBagchi et al.l ()201lD for comparisons, as do they. In 
Fig. [To] and Table U we show and calculate the correla- 
tions between our values for F and their MSP population 
results. Our statistical correlation tests indicate a very 
strong correlation between F and the number of MSPs in 
a GC, with null- hypothesis probabilities of no correlation 
below 0.013. 

5.3. j-ray fluxes 

The Fermi 7-ray Space Telescope's Large Area Tele- 
scope's unprecedented sensitivity and spatial resolution 
to GeV 7-rays have allowe d detection of num erous ra- 
dio MSPs as 7-ray sources (|Abdo et al.ll2009clia[ ). show- 
ing charac teristic hard GeV s pectra with cutoffs around 
1-3 GeV (|Abdo et al.l l2009aD . Fermi has recently de- 
tected gamma-ray emission from se veral glo bular clus- 
ters, including 47 Tuc and Terzan 5 (|Abdo et al.,,2009be 



iKong et al.l 120101: lAbdo et al.l 120101 ) . showing similar 7- 
ray spectra as radio MSPs, indicating that the observed 
7-ray flux is likely due to a population of 7-ray-emitting 
MSPs. In many clusters, no periodicities have been iden- 
tified in the 7-ray emission, indicating that numerous 
MSPs contribute to the total emission, and thus that 
measurements of the total 7-ray flux can be used to esti- 
mate the number of MSPs in the cluster. However, NGC 
6624 shows a counter- example, where a single MSP dom- 
inates the 7-ray fiux (jFreire et all 1201 J ), indicating that 
this method of estimating MSP numbers has limitations. 
Recent claims of detections of 7-ray fluxes from globu- 
lar clusters have been made for 7-ray sources lying well 
outside the half-mass radius of clusters, at low signifi- 
cance, and without evi denc e of spectral similarity to ra- 
dio MSPs ( Tam ct al. . 20TTI ). We do not trust that these 
7-ray sources represent the MSP poptflation of these GCs 
and therefore choose to evaluate the effects of our calcu- 
lations of F on the corre lations between int egrated 7-ray 
fl ux and F dis c ussed by lAbdo et all (|2010i ) . 

lAbdo et"aLl (j2010f ) measured 7-ray luminosities and 
calculated F for 8 GCs to investigate for a correlation. 
Using their reported values for 7-ray luminosities and 
our estimates for F, we find evidence (i.e., the prob- 
ability that such a correlation occurs randomly is less 
than 10 %) for a correlation between the two parameters 
fFig. [TT| Table m , in agreement with the conclusions of 
lAbdo et all ((2OIOI) . 

6. CONCLUSIONS 

In this paper we calculated the stellar interaction rate 
F for Galactic globular clusters, directly dcprojccting ob- 
served surface brightness profiles and then calculating 
F oc / p^/cr. Previous calculations have used simplified 
relations such as Fi cx p^r^/a, F2 oc pl'^r^, or have as- 
sumed King-model structures to perform integrations of 
F for clusters. Although our results are generally sim- 
ilar to previous analyses, we find significant differences 
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NGC6440 



34.0 
10' 



10^ 10^ 
stellar encounter rate 



Figure 11. Gamma-ray flux measurements from Abdo et al. 2010 
(7) versus our calculated F values (F). Core-collapsed clusters are 
shown in red. A correlation is suggested, and tabulated in Table 
131 The normalization of F is scaled so that Fngc 6121 = 1- 



Parameter 


XRBs^ 


Recycled PSs^ 


7-ray flux'' 


Pearson r 


0.942 


0.745 


0.589 


p(> |r|) 


4.5 X 10"^ 


0.013 


0.124 


Spearman r 


0.770 


0.863 


0.670 


p(> |r|) 


0.009 


0.001 


0.068 


Kendall r 


0.600 


0.674 


0.588 


P(> \r\) 


0.016 


0.006 


0.059 



Table 4 

Results of statistical tests for correlations of several different 
measurements of close binaries with our calculations of F. For all 
of these tests, the p values show the probability that a correlation 
arises randomly. Given such low prob abilities, ther e is cle ar 
evidence of correlations in all cas es. 1: IPoolev et al.l II2003I ). 2: 
IBagchi et all IMTlT) . 3; lAbdo et al.N20ia) 



in several cases, particularly for core-collapsed clusters, 
which we treat for the first time in the same way as non- 
core-collapsed clusters. A major advance in this work is 
the calculation of uncertainties in our final T estimates, 
by using Monte-Carlo simulations to incorporate the ef- 
fects of observational uncertainties. 

Comparing our F calculations with observations of 
close binaries produced by stellar interactions, we found 
strong evidence for correlations. This is in agreement 
with most previous work, but we do find significant dif- 
ferences with key recent resu lts. Comparing our F to the 
numb ers of XRBs in a GC (|Poolev et aLll^Ol iFregeaul 
I2008D . there is a suggestion that core-collapsed clusters 
may have fewer XRBs th an other GCs of similar F, in 
contrast to IFregeaul (|2008D . Comparing F to the number 
of MSPs in GCs, we find extrem ely strong correlations, 
in contrast to IBagchi et al1(|20Tl[) . Finally, we found evi- 
dence for correlation of F with the to t al 7-r ay fluxes from 
GCs, in agreement with lAbdo et al.l ()2010[) . 
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Name 




-S 


+5 


2 — i 


-<5 




+<5 




-S 




+<5 


Terzan 5 


6.80E+3 


3.02E+3 


1.04E+3 


1.86E+3 


9.34EH 


-2 


1.99E+3 


1.40E+3 


2.85EH 


-2 


3.23E+2 


NGC 7078 


4.51E+3 


9.86E+2 


1.36E+3 


5.01E+3 


2.77EH 


-2 


3.00E+2 


6.46E+3 


1.76EH 


-2 


1.66E+2 


NGC 6715 


2.52E+3 


2.74E+2 


2.26E+2 


2.55E+3 


1.33EH 


-2 


1.05E+2 


2.03E+3 


4.35EH 


-1 


5.11E+1 


Terzan 6 


2.47E+3 


1.72E+3 


5.07E+3 


1.78E+3 


8.51EH 


-2 


2.34E+3 


1.30E+3 


2.41EH 


-2 


2.79E+2 


NGC 6441 


2.30E+3 


6.35E+2 


9.74E+2 


2.56E+3 


1.84EH 


-2 


1.80E+2 


3.15E+3 


1.03EH 


-2 


1.08E+2 


NGC 6266 


1.67E+3 


5.69E+2 


7.09E+2 


2.02E+3 


1.86EH 


-2 


1.71E+2 


2.47E+3 


8.15EH 


-1 


7.78E+1 


NGC 1851 


1.53E+3 


1.86E+2 


1.98E+2 


1.54E+3 


7.66EH 


-1 


6.98E+1 


1.91E+3 


5.41EH 


-1 


5.48E+1 


NGC 6440 


1.40E+3 


4.77E+2 


6.28E+2 


1.54E+3 


5.43EH 


-2 


1.02E+3 


1.75E+3 


1.47EH 


-2 


1.55E+2 


NGC 6624 


1.15E+3 


1.78E+2 


1.13E+2 


1.20E+3 


1.13EH 


-2 


1.51E+2 


1.08E+3 


2.61EH 


-1 


2.45E+1 


NGC 6681 


1.04E+3 


1.92E+2 


2.67E+2 


9.81E+2 


6.04EH 


-1 


6.59E+1 


9.64E+2 


2.75EH 


-1 


2.63E+1 


NGC 104 


l.OOE+3 


1.34E+2 


1.54E+2 


l.OOE+3 


4.81EH 


-1 


4.64E+1 


l.OOE+3 


2.85EH 


-1 


3.08E+1 


NGC 5824 


9.84E+2 


1.55E+2 


1.71E+2 


9.16E+2 


4.15EH 


-1 


4.89E+1 


1.22E+3 


3.18EH 


-1 


2.98E+1 


Pal 2 


9.29E+2 


5.55E+2 


8.36E+2 


1.18E+3 


4.43EH 


-2 


8.02E+2 


4.57E+2 


4.89EH 


-1 


4.27E+1 


NGC 2808 


9.23E+2 


8.27E+1 


6.72E+1 


1.15E+3 


9.79EH 


-1 


l.llE+2 


1.21E+3 


2.49EH 


-1 


2.67E+1 


NGC 6388 


8.99E+2 


2.13E+2 


2.38E+2 


9.53E+2 


7.52EH 


-1 


7.41E+1 


1.77E+3 


4.89EH 


-1 


4.41E+1 


NGC 6293 


8.47E+2 


2.39E+2 


3.77E+2 


9.18E+2 


1.51EH 


-2 


2.24E+2 


1.22E+3 


3.26EH 


-1 


3.00E+1 


NGC 362 


7.35E+2 


1.17E+2 


1.37E+2 


8.09E+2 


3.61EH 


-1 


3.31E+1 


5.69E+2 


1.49EH 


-1 


1.56E+1 


NGC 6652 


7.00E+2 


1.89E+2 


2.92E+2 


7.03E+2 


1.29EH 


-2 


3.60E+2 


8.05E+2 


2.25EH 


-1 


2.13E+1 


NGC 6284 


6.66E+2 


1.05E+2 


1.22E+2 


7.50E+2 


9.82EH 


-1 


1.22E+2 


7.97E+2 


1.95EH 


-1 


1.76E+1 


NGC 6626 


6.48E+2 


9.11E+1 


8.38E+1 


6.43E+2 


1.03EH 


-2 


1.28E+2 


6.88E+2 


1.92EH 


-1 


1.79E+1 



Table 5 

F calculations and 1-a error estimations based on different equations, all normalized assuming rjvGCl04 = 1000. A portion is shown here, 

the complete table is available in the online edition of the journal. 



